Statistical mechanics of Floquet systems with regular and chaotic states 
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We investigate the asymptotic state of time-periodic quantum systems with regular and chaotic 
Floquet states weakly coupled to a heat bath. The asymptotic occupation probabilities of these 
two types of states follow fundamentally different distributions. Among regular states the proba- 
bility decreases from the state in the center of a regular island to the outermost state by orders of 
magnitude, while chaotic states have almost equal probabilities. We derive an analytical expression 
for the occupations of regular states of kicked systems, which depends on the winding numbers of 
the regular tori and the parameters temperature and driving frequency. For a constant winding 
number within a regular island it simplifies to Boltzmann-like weights exp(— /3eff^^^), similar to 
time-independent systems. For this we introduce the regular energies ^^^ of the quantizing tori 
and an effective winding- number-dependent temperature l//3eff, different from the actual bath tem- 
perature. Furthermore, the occupations of other typical Floquet states in a mixed phase space are 
studied, i.e. regular states on nonlinear resonances, beach states, and hierarchical states, giving rise 
to distinct features in the occupation distribution. Avoided crossings involving a regular state lead 
to drastic consequences for the entire set of occupations. We introduce a simplified rate model whose 
analytical solutions describe the occupations quite accurately. 

PACS numbers: 05.45.Mt, 05.30.-d, 05.70.Ln 



I. INTRODUCTION 



The response of a dynamical system to a time-periodic 
driving force is ubiquitous in both classical and quantum 
mechanics and plays a fundamental role in many phys- 
ical and technical applications. It opened the field for 
the coherent control of atoms and molecules [1], the op- 
timal control of chemical reactions [2| , or the manipula- 
tion of semiconductor-nanodevices and heterostructures 
in solids [3]. Under realistic, nonidealized conditions real 
physical systems interact with their environment. If the 
environment contains a vast number of degrees of free- 
dom, the full dynamics of the composite system is not 
traceable. The system is then interpreted as an open 
subsystem in mutual contact with a heat bath and the 
dynamics of the subsystem is characterized by its reduced 
density operator p. To evaluate the evolution of p for an 
open quantum system in a time-varying, strong exter- 
nal field is a nontrivial task, as it is permanently driven 
out of equilibrium. For only very few systems exact an- 
alytical solutions of the damped dynamics are feasible, 
in particular a driven two-level system [4] and a driven 
harmonic oscillator [o , ^ . In general systems it is studied 
numerically, e.g. with focus on tunneling, see Ref. [7] and 
references therein. Especially in the regime of weak in- 
teraction with the environment, standard methods, origi- 
nally estabhshed for time-independent quantum systems, 
have been adapted to the demands of time-periodic sys- 
tems ^^. 

The final state of the relaxation process has so far not 
received as much attention as transient phenomena, al- 
though this can be ranked as even more fundamental 
and is in fact a core question of statistical mechanics. 
The usual thermodynamic concepts for the equilibrium 
state of time-independent systems are not appHcable, 



such as the canonical distribution of Boltzmann weights, 
reached in the stationary limit of a time-independent sys- 
tem that is weakly coupled to a heat bath. The Boltz- 
mann weights e~^^^ of the eigenstates are unique func- 
tions of the eigenenergy with the temperature 1//3 = ksT 
of the heat bath as the only relevant parameter, whereas 
microscopic details of the weak coupling play no role. 
Such a stationary limit, in the sense of convergence to 
time-independent values for all dynamical variables, is 
not encountered in a periodically driven system, where 
energy is permanently exchanged between the driven sys- 
tem and the environment. Instead, the relaxation process 
finally leads to an asymptotic state that adopts the peri- 
odicity of the driving, and that in general depends on the 
microscopic details of the coupling. The density opera- 
tor of the time-periodic subsystem is best represented in 
the Floquet state basis. The Floquet states are quasi- 
periodic solutions of the Schrodinger equation for the 
time-periodic Hamiltonian without the coupling to the 
environment. In the Floquet basis the evolution equation 
for the density matrix can be approximated within the 
Floquet-Markov approach ^SrtJJ] by a Markovian quan- 
tum master equation. In the long-time limit of the evo- 
lution the Floquet states are populated with asymptotic 
occupation probabilities that can be determined from a 
system of rate equations. Beyond the numerical evalua- 
tion of such a master equation, an intuitive understand- 
ing of these Floquet occupations is still lacking. 

A related problem occurs at avoided crossings, which 
are ubiquitous in the quasienergy spectra of generic Flo- 
quet systems, since the quasienergies are bounded within 
a finite interval, < £ < ficj, where uo = 27r/r is the 
driving frequency and r is the driving period. As a con- 
sequence, the number of avoided crossings grows without 
limit for increasing Hilbert-space dimension, leading to 
a breakdown of the adiabatic theorem \13\ . In the pres- 



ence of a heat bath this problem is approached in Ref . [l2| 
where it is shown that the reduced density operator p is 
not affected by a small avoided crossing, provided that it 
is smaller than a specific effective coupling parameter and 
so is not 'resolved' by the heat bath. These findings jus- 
tify the unevitable truncation of the, in general, infinite 
Hilbert space dimension in numerical implementations. 

One way to tackle the general challenge of finding the 
Floquet occupations beyond their numerical evaluation 
is to study the semiclassical regime of one-dimensional 
driven systems. In their classical limit regular and 
chaotic motions generically coexist. This is most clearly 
refiected in phase space, where regular trajectories evolve 
on invariant tori and chaotic trajectories fill the remain- 
ing phase-space regions. According to the semiclassical 
eigenfunction hypothesis [14] almost all Floquet states 
can be classified as either regular or chaotic, provided 
that the phase-space regions are larger than the Planck 
constant. The regular states localize on the regular 
tori and the chaotic states typically spread out over the 
chaotic region. For a driven particle in a box coupled 
to a heat bath the Floquet occupations of regular and 
chaotic states were found to follow different statistical 
distributions [10]. The regular states, which in this ex- 
ample differ only slightly from the eigenstates of the un- 
driven system, carry almost Boltzmann weights, whereas 
all chaotic states have nearly the same occupation prob- 
ability. 

In this paper we concentrate on situations characteris- 
tic for strong driving, where both phase space and Flo- 
quet states are strongly perturbed compared to the orig- 
inally time-independent system. We demonstrate that 
the Floquet occupations of the states in a regular is- 
land under these conditions deviate considerably from 
the Boltzmann result. For kicked systems, making use 
of some reasonable assumptions, we derive an analytical 
expression for the regular occupations. In many cases it 
can be well approximated by weights of the Boltzmann 
type e~^^ff^^^ . This requires the introduction of the reg- 
ular energies ^^^, which are semiclassical invariants of 
the quantizing tori of the regular island, and the parame- 
ter l//3eff, which is an effective temperature depending on 
the winding number of the regular island. Furthermore, 
we give an overview and interpretation for the occupa- 
tions of other typical Floquet states in a mixed phase 
space, such as states on a resonance island chain, beach 
states, and hierarchical states. 

Avoided crossings in the Floquet spectrum can lead to 
severe changes in the occupations if they are larger than 
an effective coupling parameter. The effect can be intu- 
itively explained by a set of effective rate equations with 
an additional rate R^'^ between the states of the avoided 
crossing [12]. It can be exploited for a switching mech- 
anism in driven quantum systems, e.g. a weakly driven 
bistable system [15]. For the above occupation distribu- 
tions for regular and chaotic states we demonstrate dras- 
tic consequences if a regular state has an avoided crossing 
with either a regular or a chaotic state. We introduce a 



simplified rate model whose analytical solution describes 
the Floquet occupations accurately. 

The paper is organized as follows: in Sec. [IT] the mi- 
croscopic model of driven dissipative systems and the 
Floquet-Markov description of its asymptotic state are 
sketched and the relevant coupling operators are intro- 
duced. Section |TTl] presents general occupation charac- 
teristics for the example of a driven quartic oscillator 
(Sec. nil Ap and a kicked rotor (Sec. IIIIBp . They are re- 
lated to the corresponding rate matrices (Sec. IIIIC]) . In 
SecHVlwe derive an analytical expression for the regular 
occupations of kicked systems, depending on the wind- 
ing numbers of the regular tori. For a constant winding 
number a simplification to the Boltzmann-like weights 
g-/3eff£;^ is shown (Sec. lIVCp . An example where this 
is not possible is also discussed (Sec. lIVPp . Section IVl 
gives an overview of the occupation characteristics of 
other types of Floquet states. In Sec.|Vl]the infiuence of 
avoided crossings on the Floquet occupations is demon- 
strated, which are compared to the analytical solutions 
of a simplified rate model (Sec. I VI Ap . Section |VlI] sum- 
marizes the paper. 



II. MASTER EQUATION IN TIME-PERIODIC 

SYSTEMS 



The coupling of a quantum system with the Hamilto- 
nian Hs(t) to a heat bath is modelled in a standard way 
by the composite Hamiltonian [ 16l] 



H{t) = Hs{t)^H,^Hsb. 



(1) 



the 



bath Hamiltonian Hi^ = 

"^ ' describes an ensemble of noninter- 



Herein, 

T (^ 

acting harmonic oscillators coupled via the interaction 
Hamiltonian Hsh to the system. In spatially extended 
systems this interaction is commonly assumed to be 
bilinear. 



J^sh — / V ^ri'^ni 



(2) 



with some coupling operator A of the system. The 
properties of the system-bath coupling are speci- 
fied by the spectral density of the bath J{uo) := 

2 

f Z)n t;^^ [^(^ ~ ^n) - K^ + ^n)]. In the continuum 
limit the spectral density is assumed to be a smooth 
function which is linear for an Ohmic bath. An ex- 
ponential cutoff beyond the spectral mode uoc leads to 
J(uj) = 7^co'e~l^l/^^, where rj is proportional to the clas- 
sical damping coefficient. 

In the absence of the heat bath the solutions of the 
time-dependent Schrodinger equation for the isolated sys- 
tem with the r-periodic Hamiltonian 



H,,{t + T)=H,{t) 



(3) 



are the Floquet states \ilji{t)). These can be factorized 
into a product 



mt)) 



-isit/h\ 



Ui(t)) 



(4) 



of a phase factor with the quasienergy £j and a periodic 
state vector \ui{t)), 



\Ui{t + T)) = \Ui{t)), 



(5) 



with the period r of the Hamiltonian. 

In the presence of the heat bath the state of the sys- 
tem is described by the reduced density operator p(t). 
Its equation of motion for time-periodic quantum sys- 
tems has been derived within the Floquet- Markov ap- 
proach |8Ml2]: herein the Floquet formalism ensures a 
non-perturbative treatment of the coherent dynamics of 
the driven system. The density operator is represented 
in the set of the time-periodic state vectors \ui{t))^ which 
form a complete ort honor mal basis at all times t. The 
coupling to the heat bath is treated perturbatively in 
second order of Hgh^ which is valid in the limit of weak 
coupling between the driven system and the bath. This 
approximation requires a rapid decay of bath correla- 
tions compared to the typical relaxation time of the sys- 
tem. In this paper we use a cutoff frequency uoc — lOOcj, 
which is large compared to the frequency uj = 27r/r 
of the driving. In the following we restrict the discus- 
sion to the limit of large times, much larger than the 
relaxation time. In this limit the density-matrix ele- 
ments pij = (ui(t)\p(t)\u^(t)) are approximated as time- 
independent [qI, Ii2| • Note that the corresponding den- 
sity operator, ^^ • \ui{t))pij{uj{t)\^ is still time-periodic 
because of the inherent time-dependence of the \ui{t)). 
In this paper we restrict to the weak-coupling regime, 
where the system-bath coupling is small compared to all 
quasienergy spacings of a truncated Hilbert space (see 
discussion below). The Floquet occupations pi = pu then 
obey the set of rate equations 



= -PiY^ Rij + ^PjRji, 



(6) 



The correlation function g{E) = 7r~^ n 13(E) J (E/h) of the 
bath coupling operator contains the spectral density J{uj) 
and the thermal occupation number n/3 (E) of the boson 
bath with temperature 1//3. 

The reduction to the set of rate equations (J6]) seems 
possible only for systems with a finite dimension of the 
Hilbert space, since otherwise the quasienergies densely 
fill the interval [O^hu). However, as demonstrated in 
Ref. [12|, near degeneracies much smaller than the cou- 
pling strength are not resolved by the interaction to the 
heat bath and do not infiuence the asymptotic density 
operator. The Hilbert dimension can therefore be trun- 
cated, keeping only those Floquet states of non- negligible 
occupation. 

Equation (J6]) is formally identical to the familiar sys- 
tem of rate equations describing the equilibrium state 
in time-independent systems. In contrast, however, 
specifics of the time-periodic system are present in the 
rates, Eq. (O, whose structure does in general not allow 
a detailed balance relation |11|. 



A. Coupling operator for extended and for cyclic 
systems 

For extended systems we assume as usual [l6| the linear 
coupling operator 



A = X 



(10) 



in Eq. (|2j) for the interaction Hamiltonian Hsb with the 
heat bath. 

For cyclic systems defined on the unit interval [0,1) 
with periodic boundary conditions in x the coupling op- 
erator X would be discontinuous at the borders of the 
interval and the coupling to the heat bath would there- 
fore not be homogeneous. An adapted coupling scheme 
with the interaction Hamiltonian 

Hsb = / da — sm{27rx -\- a) y^ CnXnS {a — an) (11) 
^0 27r ^ 



which are independent of the damping coefficient 77. Note 
that the rate equations beyond this weak-coupling regime 
would also contain the non-diagonal elements pij {i ^ j). 
The rates 



K 



that describe bath-induced transitions between the Flo- 
quet states, use the Fourier coefficients 



'^ Jo 



of the time-periodic matrix elements 



(8) 



(9) 



has been proposed in Ref. [17] for such situations. The 
angles an characterize the individual bath oscillators and 
are equidistributed in the interval [0, 27r). The new spec- 
tral density 



^^— (5(a-an) [S{UJ - UJn) - S{UJ +CJn)] 



Jiuj^a) = — 7 

2 ^ rrinUJn 

(12) 
hence factorizes into independent factors, J{uj^a) = 
J(co')/(27r), with the spectral density J{uj) as defined 
above and the homogeneous angular density l/(27r). The 
spatially periodic interaction Hamiltonian (pT]) is con- 
tinuous and, by virtue of the equidistributed angles 
a^, models the interaction with a homogeneous envi- 
ronment, where no position is singled out. Making use 
of the trigonometric addition theorem sin(27rx -\- a) = 



sin(27rx) cos(Qf) + cos(27rx) sin(Qf), the interaction Hamil- 
tonian (pT]) leads to the same system of rate equations, 
Eq. (J6]). But now the rates 



Ri 



d(1) , d(2) 



(13) 



are composed of two independent contributions from the 
simpler coupling operators 



A^^^ = sin(27rx)/(27r), 
A^^^ = cos(27rx)/(27r), 



(14) 
(15) 



respectively, to be used in the interaction Hamiltonian 
Hsi) in Eq. ([2]). For this result we use that the mixed 
second-order term J^ ^ da J^^ da' cos(Qf) sm{a')6{a — 

an)S{a' — an) = J^^ da cos(a) sin(a) = vanishes, while 
the other two terms give rise to the operators in Eqs. ([H]) 
and ([T5l). 



OCCUPATIONS OF REGULAR AND 
CHAOTIC FLOQUET STATES 



III. 



the driving period, r = 27r/co', features a chaotic domain, 
see Fig. [IJa). Furthermore, there are two distinct regu- 
lar regions: first the highly excited tori, which are only 
slightly influenced by the driving, and second a regular 
island embedded in the chaotic sea. 

The Floquet states are determined using the it^t')- 
technique [18]. Their energy expectation value, 

{^lji{t)\Hs{t)\^lJi{t)) = {ui{t)\Hs{t)\ui{t)), oscillates with 
the period of the driving. It is convenient to introduce 
the cycle- averaged energy 

(E^) :=- [ \t'{u,{t')\Hs{t')\u,{t'))-E^. (18) 
^ Jt 

An energy shift Eq is determined by the classical periodic 
orbit at the center of the regular island, such that there 
the cycle- averaged energy is zero. 

For a sufficiently small value of h almost all Floquet 
states can be classified as either regular or chaotic ac- 
cording to the semiclassical eigenfunction hypothesis [14]. 
The regular states are localized on the regular island 
and can be ordered by a quantum number m, whereas 



In this section we study the Floquet occupations for 
two classes of periodically driven systems, the additively 
driven quartic oscillator as a representative of a continu- 
ously driven system and the quantum kicked rotor from 
the class of kicked systems. In both cases we demonstrate 
that the Floquet occupations of regular and chaotic Flo- 
quet states follow tremendously different distributions. 
Similar observations have been made for a driven quan- 
tum particle in a box [10]. While there it was found 
that regular states carry Boltzmann weights, we find typ- 
ically significant deviations. Numerical results will be 
presented in this section and the quantitative analysis of 
the occupations of the regular states will be deferred to 
Sec.lIVl 



A. Driven Quartic Oscillator 

We consider as an example of a continuously driven 
system the additively driven quartic oscillator with the 
Hamiltonian 



1.0 



Hs{t) = ^ + ^0 (^ + /^— cos(a;t) 
2m V Xq xq 



(16) 



We introduce the dimensionless quantities x = x/xq, 
H = H/Vo, p = p/po with po = VrnVo, i=t- po/{xom), 
uj = uj ■ {xom)/po^ and also h = h/{xopo)^ the ratio of the 
Planck constant to a typical phase-space area. In the fol- 
lowing we omit the overtilde and then the dimensionless 
Hamiltonian reads 



Hs{t) 



2 



+ X'^ + KXCOs{u)t). 



(17) 



At Ai: = 0.2 and uj = 5/6 the stroboscopic Poincare section 
of the phase space (x,p) at integer multiples t = nr of 
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FIG. 1. (Color online) (a) Stroboscopic Poincare section 
of the classical phase space of the driven osillator, Eq. ()17p . 
The size of the chosen dimensionless Planck constant h is 
indicated in the lower right corner, (b) Floquet occupations pi 
vs. cycle- averaged energies (Ei) compared to the Boltzmann- 
like prediction exp(— /3(^i)) (dashed line). The insets show 
Husimi representations of a regular Floquet state localized in 
the central island, a chaotic Floquet state, and a regular state 
on a surrounding torus. The parameters are k = 0.2, u = 5/6, 
h = 0.002, and /3 = 100. 



the chaotic states typicahy spread over the entire chaotic 
phase-space area and fluctuate irregularly. The different 
types of states are visualized in phase space by means 
of their Husimi representation H^{x^p) := \{a{x^p) \^p)\, 
i.e. the projection onto the coherent states \a{x^p)) cen- 
tered at the phase space points (x,p), see insets in 
Fig. [IJb). In the example of Fig. [1] the small central 
island is about 26 times larger than the dimensionless 
Planck constant /i, indicated in the lower right corner of 
Fig. [TJa). Thus the central island supports 26 regular 
states. Besides, there are 97 chaotic states, spreading 
over the chaotic region. It is surrounded by regular tori, 
on which a further group of infinitely many regular states 
are localized. The cycle- averaged energies (Ei) of these 
three different types of Floquet states form distinct inter- 
vals with only small overlaps, as indicated by the arrows 
in Fig. [IJb). The regular states of the central island are 
lowest in cycle- averaged energy, followed by the chaotic 
states and finally the high excited regular states of the 
surrounding tori. 

The oscillator coupled to a heat bath is treated as ex- 
plained in Sec.|II]with the linear coupling operator A = x 
for this spatially extended system. We restrict the con- 
sideration to the weak-coupling regime, where the occu- 
pations are well described by the rate equation ([6]). In 
Fig. [IJb) the resulting Floquet occupations pi are shown 
as functions of the cycle- averaged energies (Ei). The 
monotonously falling occupations at low values of (Ei) 
belong to the central regular island, with the state in the 
center of the island having the highest occupation. At 
intermediate values of (Ei) one finds the occupations of 
the chaotic states, that fluctuate around a mean value 
Pch, with a very small variance compared to the range 
of occupations of the regular states. This is similar to 
the observation for chaotic states in Ref. [10] for a driven 
particle in a box. At high values of (£^i), there are again 
monotonously falling occupations belonging to the regu- 
lar states of the surrounding tori. 

The observed characteristics of the occupations pi are 



0.5 



clearly in contrast to the naive expectation pi ~ e 



motivated by the Boltzmann weights of equilibrium ther- 
modynamics. Note that even the occupations of the (low- 
energy) regular states notably differ from the Boltzmann 
result, indicated by the dashed line in Fig.fTJb). A quan- 
titative analysis of these observations will be presented in 
Sec. II VI for the numerically more convenient kicked rotor. 



B. Kicked rotor 

Kicked quantum systems feature all essential phase- 
space characteristics of periodically driven systems. They 
allow for a simplified numerical and conceptual treat- 
ment. As a paradigmatic model for a driven system with 
a mixed phase space we consider here the quantum kicked 
rotor. Its dynamics is generated by the Hamiltonian 

Hs{t) = T{p) + V{x) • r ^ J(t - nr), (19) 
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FIG. 2. (Color online) (a) Stroboscopic Poincare section of 
the classical phase space of the kicked rotor, Eq. ([20]) . for 
K = 2.9. The size of the chosen dimensionless Planck con- 
stant h is indicated in the lower right corner, (b) Floquet oc- 
cupations pi vs. cycle- averaged energies (Ei) compared to the 
Boltzmann- like prediction exp(—/3(£^i)) (dashed line). The in- 
sets show Husimi representations of two regular and a chaotic 
Floquet state, (c) Floquet occupations pm of regular states 
vs. regular energies E^^, defined in Eq. ([2T]) . compared to the 
Boltzmann-like prediction Eq. (|40p with the inverse effective 
temperature /3eff = 0.85/3 of Eq. (|39p (red solid line), com- 
pared to the inverse bath temperature /3 (dashed line). The 
parameters are h — 1/210 and /3 = 100. 



with the kinetic energy T{p) = p^ /{2m) and the potential 
V{x) = Fo cos(27rx/xo) acting at the kicks. We study it 
on a two- tor us with dimensions xq and po = vfix^jr. We 
make the Hamiltonian dimensionless by a similar trans- 
formation as in Sec. IIII A| where now H = H ■ m/pQ and 
the dimensionless kick period is r = 1. With the rescaled 



kick strength k. = Vq - (27r)^m/pQ the Hamiltonian reads 



Hs{t) 



{2ny 



■ cos{27rx) ^ 5 {t-n). (20) 



In an intermediate regime of the kick strength the 
Poincare section of the phase space features a regular 
island embedded in the chaotic sea, see Fig. [Sfa) for 
K. = 2.9. 

The Floquet states are evaluated as eigenstates of the 
time evolution operator U over one period r, which 
factorizes into a potential and a kinetic part U = 
^-irVix)/h^-irT{p)/h^ r^^^ quantization on the two-torus 
relates the effective Planck constant h to the dimension 
N of the Hilbert space by the condition h = 27Th = 1/N. 
For /i = 1/210 the area of the regular island supports 23 
regular states. 

The asymptotic state of the kicked rotor weakly cou- 
pled to a heat bath is again determined from Eq. (J6j), 




210 



?(2) 



with the composite rates Rij = i?- • +i?- • from Eq. ([13]) FIG. 3. Density representation of the rate matrix Rij, with 



appropriate for a cyclic system. The resulting Floquet 
occupations pi are shown in Fig. [2fb) as functions of the 
cycle- averaged energies (£^i), with Eq = — n / (27:)'^ . The 
regular and chaotic states again are ordered with respect 
to this quantity, as indicated by the arrows in Fig. [2fb). 
The regular states have small values of {E)^ since both 
kinetic and potential energies are minimal in the center 
of the regular island, whereas the chaotic states have a 
stronger overlap with regions of phase space with higher 
energies. Similarly as for the driven oscillator, the regu- 
lar occupations depend monotonously on (£^), while the 
occupations of the chaotic states seem uncorrelated with 
the cycle- averaged energy {E) and form a plateau with 
only weak fluctuations around a mean value p^h- 



C. Rate matrix 

The Floquet rate matrix Rij determines the Floquet 
occupations via Eq. ([6]). In Fig. [3] we show Rij for both 
the driven oscillator in Fig. [1] (Fig. [3fa)) and the kicked 
rotor in Fig. [2](Fig. [3l^b)). Employing {Ei) as the order- 
ing parameter for the entries z, j, the regular and chaotic 
parts are well-separated, revealing a distinct block struc- 
ture of the matrix [lO|. There are only few rates be- 
tween the regular and the chaotic subspaces. Similarly, 
the rates between the two different regular subspaces in 
the case of the driven oscillator in Fig. [3fa) are practi- 
cally zero. Also by virtue of the chosen ordering, the 
regular domains feature a band structure with particu- 
lar dominance of the first off diagonals (nearest-neighbor 
rates). The rates in the subspace of the chaotic states on 
the contrary fluctuate strongly. 

One can thus observe a close relation between the 
structure of the rate matrix Rij and the resulting set 
of occupations. First, the almost independent behaviour 
of the occupations of regular states and chaotic states 



entries i ^ j sorted by increasing {Ei) for (a) the driven os- 
cillator, Eq. ([T7|, and (b) the kicked rotor, Eq. ([20)), with 
enlarged domain of regular island states, revealing strong 
nearest-neighbor rates. For parameters see Figs. [T] and [21 
respectively. 



owes to the relatively weak rates Rij connecting the cor- 
responding subspaces. Furthermore, the random charac- 
ter of the chaotic rate submatrix gives rise to the equally 
random character of the set of chaotic Floquet occupa- 
tions. Note that for the kicked rotor we find in the semi- 
classical limit /i ^ 0, that the mean value Pch decreases, 
since more and more regular states emerge. Also the rel- 
ative variance {pi — pch) /Pch ^f the chaotic occupations 
Pi decreases in this limit and we observe a universal scal- 
ing in the semiclassical limit which can be analyzed with 
the help of a random-rate model [19|. These aspects of 
the chaotic occupations are not explored in this paper, 
instead the focus is set on the regular occupations. 



IV. REGULAR STATES 

The observations in Figs. [1] and [2] indicate that the 
asymptotic state of a time-periodic system in weak inter- 
action with a heat bath carries signatures of the classical 
phase-space structure. The Floquet occupations of the 
regular and the chaotic states behave very differently, 
e.g. as functions of the cycle-averaged energy {E). In 
this section we focus on the asymptotic occupations of 
the regular states, which we label by their quantum num- 
ber m starting with tti = for the state in the center of 
the island. Figures [TJb) and [2fb) suggest a roughly ex- 
ponential dependence for the regular occupations Pm as 
functions of the cycle- averaged energies {Em)- However, 
the regular occupations are different from the Boltzmann 



weights e~^^^"^^ with the true inverse bath temperature 
/3. In fact, there is no physical reason for a coincidence 
with the Boltzmann distribution when expressed in terms 
of the, quahtatively suitable but arbitrary, energy mea- 
sure (Em) [11]. In the following sections we therefore 
make use of an alternative energy measure for the regu- 
lar states, the regular energy E^^ fSec. IIV A|) , allowing 
us to consistently parametrize the regular occupations as 
functions of E^^ (Sec. lIVBp . Often, this functional de- 
pendence is approximately exponential (Sec. IIVU]) . Ex- 
amples are presented in Sec. IIVDI 



A. Regular energy E^^ 

A time-periodic system is equivalent to an autonomous 
system with the time as an additional coordinate, leading 
to the Hamiltonian H'^ix^p^t^pt) = E[s{x^P'-,t) -^ Pt in 
the extended phase space, which has periodic boundary 
conditions in t. This allows the application of Einstein- 
Brillouin-Keller(EBK)-quantization rules for the regular 
tori and the determination of semiclassical Floquet states 
on the quanti zing tori and their associated semiclassical 
quasienergies |2Q|, |21|. 

We introduce the regular energies 



Ei:^ := huji 



m 



W™ + w. 



(21) 



Herein, Vrn is the winding number, i.e. the ratio of the 
winding frequency of a trajectory on the mth torus 
around the central orbit to the driving frequency uo. 
Furthermore, {L)^ is the long-time average of the La- 
grangian L = px — Hg for an arbitrary trajectory on the 
mth torus. For convenience we add the time-averaged 
Lagrange function {L)c of the central orbit of the island. 
The regular energies of Eq. (|2T]) are related to the semi- 
classical quasienergies [20, 21] by 



^m = E]^^ - {L)c mod huj, 



(22) 



whereas there is no relation to the cycle- averaged ener- 
gies (Em)- The time- averaged Lagrangian (L) varies only 
slowly inside the island. The winding number u likewise 
varies slowly across the island. To determine u we have 
applied the frequency map analysis [22| being based on a 
Fourier decomposition of the quasiperiodic orbits within 
a stable regular island. Note that due to nonlinear reso- 
nances and small chaotic layers within a regular island 
the semiclassical quantization might require interpola- 
tions of the quantities Urn and {L}rn or the introduction 
of a fictitious integrable system [2l|, |23| . 

Figure [2fc) shows the occupations pm of the regular 
states of the kicked rotor as functions of regular ener- 
gies E^^^. The functional dependence of the occupations 
Pm is close to exponential, but also different from the 



Boltzmann weights e 



-^Ei^: 



with the true bath temper- 




FIG. 4. (Color online) Ratio between the occupations pm 
and an exponential fit pQe~^^''^^^'^~^^^ for (a) the kicked rotor 
and (b) the continuously driven oscillator using the regular 
energies Cm = ^m^ (red circles) and the cycle- averaged ener- 
gies Cm = {Em) (black diamonds). For parameters see Figs. [2] 
and[T] respectively. 



(Em)- This is demonstrated in Fig. HI where the ratio 
between the occupations Pm and the respective exponen- 
tial fit pQe~^fi*(^"^~^o) is shown for Cm being the regular 
energy E^^ (red circles) and the cycle-averaged energy 
{Em) (black diamonds). The fit involves the parameter 
/^fit •= (logpi — logpo)/(eo — ei). For the kicked rotor, 
Fig-IHa), the considered ratio for Cm = E^^ is close to 1 
for the majority of regular states, whereas the ratio for 
Cm = {Em) systematically deviates from 1 already for 
smaller values of m. This indicates that the exponential 
scaling is far better fulfilled by using the regular energies 

E;reg 

Likewise, Fig. HJ^b) shows the same ratio for the reg- 
ular states of the central island in the driven oscillator, 
Eq. (pT|) . The regular energies are again determined ac- 
cording to the above semiclassical quantization, where 
the frequency map analysis is applied to the solutions 
of the classical equations of motion, evaluated in the 
Poincare section. Here, the quality of the fit with re- 
spect to El^^ is only marginally better as the the fit with 
respect to {Em), see Fig. SJ^b). We have evidence, that 
the existence of next-nearest-neighbor rates is responsi- 
ble for this. Note that for other examples of continuously 
driven systems we typically find a better quality of the fit 
with respect to £^^^, similar to the situation in Fig.HJ^a). 



B. Restriction to nearest-neighbor rates Rm,m±i 

In this section, the ratio of the rates Rm,m-\-i and 
Rm-\-i,m between two neighboring regular states m and 
?n + 1 is analyzed for kicked systems. With the help of 
a detailed balance condition the occupations Pm can be 
related to the winding numbers of the regular tori. 

In the lower part of Fig. [3fb) the rate matrix Rij = 



?(i) 



ature 1//3. However, the assumption of an exponential 
dependence of Pm vs. E^^ is fulfilled far better than vs. 



R\j^ + Rij , Eq. ([13]), for the regular subspace of the 
kicked rotor is shown. We remind the reader that 
the indices are ordered by increasing {E)^ coinciding 
with the natural order of growing quantum number 



m. Figure ^h) illustrates that the nearest-neighbor 
rates Rm,m±i are dominant among the regular states. 
These nearest-neighbor rates are mainly contributed by 
the rates R^, originating from the coupling operator 

^4^"^^ = sin(27rx)/(27r), whereas the rates Rm,m±2 be- 
tween next-nearest neighboring states are mainly due to 
the contribution R\-^ of the coupling operator A^^^ = 
cos(27rx)/(27r). In the following analytical considerations 

(2) 

we will neglect the contribution of R\j and in addition 
approximate the coupling operator A'^^^ = sin(27rx)/(27r) 
inside the regular island at x = 0.5 by the linear coupling 
operator A = —x. Using the resulting rate matrix in 
Eq. (J6]), we observe almost the same regular occupations 
as in Fig. [2fb), which in first approximation differ from 
the latter only by a tiny m-independent factor. Only the 
occupations of the chaotic states are strongly affected by 
the different coupling scheme, as for them the disconti- 
nuity of the coupHng operator x at the border of the unit 
cell is not negligible. 

In the rate matrix due to the linear coupling operator, 
A = —X, the nearest-neighbors dominate strongly, and 
the next-nearest-neighbor rates Rm,m±2 of regular states 
m and m ± 2 are zero for symmetry reasons. We will 
neglect higher-order rates in the following analysis. Thus, 
the total rate balance among the regular states can be 
reduced to the detailed balance condition 



Pm+l 
Pm 






(23) 



between two neighboring regular states m and m + 1. 

Assuming Eq. ([23|) and using the definition of the rates 
in Eq. ([7]) the occupation ratio between neighboring reg- 
ular states becomes 



Pm+l 
Pm 



E^\Xm{K)\^g{Cm^Khu) 



(24) 



J^xkmWI g{Cm^Khj)e 



, (25) 



where the properties A^^^+i(i^) = A'^_^-^^{—K) and 
g{E) = g{—E)e~^^ have been used, and where we intro- 
duced the short-hand notations 



:= A 



771+1, m 



(26) 



for the (regular) nearest-neighbor matrix elements of the 
operator A = —x and 



c« 



-m+l 



(27) 



in the arguments of the correlation function g{E). 

If there were just a single Fourier component Xm{K*), 
which is approximately the case for a weakly driven 
system, then the occupation ratio would simplify to 

Pm+i/Pra = e~^^^'^^^* ^^\ resulting in Boltzmann-like 



occupations. In general, however, several components K 
have to be considered. 

Multiplying numerator and denominator of the frac- 
tion in Eq. (|25|) with 1x^(0)1"^ ^ (Cm)"\ 



Pm+l 



^K 1^^(0)1^ giCm) ^ 



(28) 



we introduce ratios of the matrix elements and of the cor- 
relation functions. The ratio of the correlation functions 
reads, using their definition in Sec. HH 



g{Cm) 



1 



•exp 



\Cm\-\Cm^Khcj\ 
hWr 



(29) 



The last factor in Eq. (|29|) is close to 1 and will be omit- 
ted in the following, such that the cjc-dependence is ne- 
glected. This is possible since we are interested here in 
the case ujc ^ oj and we use the fact that only small in- 
tegers K contribute significantly to the sums in Eq. (|28]) . 
For the required ratio of the matrix elements we ap- 
proximate the evolution of the coupling matrix elements 
Xm{t) for < t < r by 



Xm{t) - Xm{t =0)e 



-iCmt/h 



1 



LKrr.r/h _ A 



(30) 

using the factorization of the time evolution opera- 
tor for kicked systems and the approximate commu- 
tation relations with the operator x on the two-torus. 



and [x, e-^^(^)^/^] ^ tr{p)e 



[x,e 

These commutation relations, which are exact in the infi- 
nite Hilbert space, apply here in very good approximation 
to the regular states, as these are almost independent 
of the periodic boundary conditions on the two-torus. 
The coupling matrix elements Xm(t) are time-periodic 
and have the Fourier components 



Xm{K) 



Xm{t = 0) 



27r2 
whose ratio simplifies to 






K 



1 — cos 27r 



\xmm' 



1 



^7; 



Finally, inserting Eqs. (|29|) and 
tion ratio of Eq. (|28|) yields 



Pm+l 
Pm 



= F 



Sm 

huj 



,phuj 



(31) 

(32) 
into the occupa- 

(33) 



with the function 



F{z,b):-- 



Ek(^ + ^) (e 



3 (JK+z)b 



^j^iK-zy^ei^-^)"-!)- 



(34) 
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It is invariant under an integer shift of the first argument, 
F{z + Ko, h) = F{z, 6), with Kq e Z [24]. We choose the 
shift Kq, such that 



-m+l 



Kohw ■■ 



EZ\i 



£.reg 



(35) 



is fulfilled, which is possible according to Eq. (|22]) . This 
allows us to replace Cm in Eq. (J33j) with the regular energy 
spacing EZ\i^m '•= ^m+i - ^m^. leading to 



Pm+l 
Pm 



/ E;reg 
77, / 7n+l,m 



, phuj 



(36) 



Based on Eq. (|2T]) we approximate this energy difference 
by the winding number 



^m+l,m 



hhJUri 



(37) 



which is exact for a harmonic oscillator-like island with 
m-independent winding number Vrn and {L)m and is a 
reasonable approximation even for more generic islands. 
The occupation ratio then becomes a function of the 
winding number. 



Pm+l 
Pm 



F {Um.Phuj) . 



(38) 



With that an analytical prediction for the occupation 
of the regular states of a kicked system is found, valid 
under the assumption of dominant nearest-neighbor rates 
Rm,m±i- It is a function of the winding number of the 
mth quantizing torus and of the parameters temperature 
and driving frequency. 



C. Assumption of constant winding number u 



The function F in Eq. ([38]) becomes independent of 
m if the winding number u is constant throughout the 
regular island. It is then appropriate to introduce an 
effective temperature l//3eff by 



/3eff:= 



logF{u,f3huj) 



(39) 




0.00 



0.25 



ly 0.50 



FIG. 5. Inverse effective temperature /3eff//3 according to 
Eq. ()39|) vs. winding number v for uJc/oj ^ 1 and three dif- 
ferent temperatures. 



With this new parameter the occupation ratios are ex- 
pressed in a form analogous to the Boltzmann weights. 



Pm+l/P7i 



-^effEZl^^^ 



(40) 



The ratio /3eff//3 is shown in Fig. [5l Its value is smaller 
than 1 and it is symmetric in v. For i^ ^ 0, where 
the kicked system approaches its static limit, the true 
bath temperature 1//3 is retained. A substantial devia- 
tion from the true bath temperature, Pef^/P "C 1, takes 
place around u ^ 0.5. 

For generic islands with non-constant winding number, 
where the regular energy spacings £^^^i ^ — huum are 
m-dependent, the exponential scaling, Eq. (|4Q|) . is still 
approximately valid if /3eff varies only moderately. Fig- 
ure [5] indicates that this is fulfilled especially good for 
values i^ < 0.2. But even beyond this interval we observe 
a good agreement of Eq. (|4Q|) with the regular occupa- 
tions. This breaks down for islands with winding num- 
bers varying close to u = 0.5, where /3eff is particularly 
sensitive to variations ofu. As /3eff develops a pronounced 
m-dependence there, the approximation (|4Q|) is no longer 
adequate and the general equation ([38|) has to be used 
instead, as demonstrated in Fig. [6] below. 



D. Examples 

Figure [2]^c) shows the regular Floquet occupations of 
the kicked rotor, Eq. ([20]) . vs. the regular energies E^^ 
and demonstrates excellent agreement with the above 
predicted exponential weights e"^^^^^^^ (red solid line) 
for almost all regular states. The regular energies E^^ 
are here slightly m-dependent, as the winding number 
decreases in the island monotonously from uq = 0.32 for 
the first (m = 0) regular state to 1^22 = 0.28 (m = 22) 
for the outermost regular state. Deviations from the 
exponential distribution occur for the outermost regu- 
lar states, m > 20, only. These have a stronger weight 
outside the regular island and are thus coupled stronger 
to the chaotic states. The non-negligible rates between 
these regular states and the chaotic states, see Fig. [3l 
enforce a gradual adaptation between the outermost reg- 
ular occupations and the occupation level of the chaotic 
states. Besides, Fig. [2fc) shows the large discrepancy of 
the Boltzmann weights e~^^^ with the true bath tem- 
perature 1//3. 

The exponential distribution, Eq. (|^Q|) . requires a con- 
stant or moderately varying winding number inside the 
island. If however v varies close to v = 0.5, where /3eff 
is particularly sensitive to variations of u, see Fig. [5l the 
exponential distribution is no longer adequate. For the 
kicked rotor this is the case for k, -^ 4, where the cen- 
tral periodic orbit bifurcates and the regular island hence 
splits into two islands. As an example. Fig. [6] shows 
the occupations for the kicked rotor for hc = 3.9, with 
uq = 0.44 and i^ig = 0.38. The regular occupations in 
Fig. ini^c) are well described by the general prediction. 
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FIG. 6. (Color online) (a)-(c) Analogous to Fig. [2] for the 
kicked rotor with tv — 3.9. The insets in (b) show Husimi rep- 
resentations of two regular states. In (c) the occupations of 
the innermost 20 regular states, which are not yet affected by 
phase- space structures at the border of the island, are well de- 
scribed by the analytical prediction (red solid line) of Eq. (|38p . 
The parameters are h — 1/400 and (5 — 500. 



tical to the Boltzmann weights e~^^^"^^ with the true 
temperature 1//3 [10,]. However, in this particular exam- 
ple the regular region is almost identical to the undriven 
system, leading to Boltzmann weights for the regular 
states by the following reasoning: the regular states of the 
driven box potential, which emerge from the highly ex- 
cited eigenstates of the undriven box, still strongly resem- 
ble the latter and change only slightly during the driving 
period r = 2t: /uo. By that, one dominant Fourier contri- 
bution i^* in the coupling matrix elements is singled out, 
Xm{K) ^ for K 7^ K^. In this situation the occupation 
ratio ([25D simplifies to Pm+i/Pm = e-^(^-+^*^^) and 
the detailed balance among the regular states, Eq. (J23j), 
is fulfilled accurately. Since at the same time their cycle- 
averaged energies are close to the eigenenergies in the un- 
driven potential, the occupations are close to the Boltz- 
mann weights. Tiny deviations for the lowest regular 
states close to the chaotic region are visible in Fig. 1 of 
Ref. [10], which we attribute to rates occurring between 
the regular and the chaotic Floquet states. 

As substantiated in this section, systematic and much 
stronger deviations from the Boltzmann behavior can oc- 
cur in generic situations, especially in situations charac- 
teristic for strong driving, where the phase-space struc- 
ture and the Floquet states are strongly perturbed com- 
pared to the original time-independent system. The 
strong driving allows us to study the Floquet occu- 
pations far from the thermodynamic equilibrium situa- 
tion, encountered in the time-independent system, and 
at the same time to have dominant regular structures 
present in the classical phase space. These host, a suf- 
ficiently semiclassical regime presumed, a series of reg- 
ular states, which under the condition of pronounced 
nearest-neighbor rates and only small rates to the sub- 
space of the chaotic states are occupied with weights 
given by Eq. (|38|) . For constant or slowly varying winding 
number u the occupations even simplify to the exponen- 
tial weights, Eq. (|4Q|) , with the effective temperature of 
Eq. ([391). 

A generic modification of the analytical predictions 
of this section, takes place as a consequence of avoided 
crossings. We will go back to this point in Sec. |Vll 



Eq. (|38|) (red Hue), but clearly deviate from the expo- 
nential approximation, Eq. (|4Q]) . 

The derivation of the analytical occupation ratios (|38|) 
and (|30|) is based on assumptions which are justified for 
kicked systems only. For continuously driven systems an 
analogous prediction remains a future challenge. Fig- 
ure [IJb) demonstrates that a significant and systematic 
deviation of the occupations from the Boltzmann result 
is observed also for continuously driven systems. 

For another continuously driven system, the driven 
particle in a box, the regular occupations are almost iden- 



V. IMPLICATIONS OF ADDITIONAL 
CLASSICAL PHASE-SPACE STRUCTURES 



The set of Floquet states in the examples of the last 
section are dominated by regular states in large regu- 
lar islands and chaotic states. Apart from these, other 
types of Floquet states can exist, depending on the struc- 
tures in the classical phase space and the size of the ef- 
fective Planck constant h. The following section gives an 
overview of the fingerprints of such additional types of 
Floquet states on the distribution of the Floquet occu- 
pations Pi. 
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A. Nonlinear resonance chains 

Apart from the islands centered at stable elliptic fixed 
points of period 1, there are nonlinear r:5-resonances con- 
sisting of r regular islands around stable periodic orbits 
of period r, see Fig. [Tfa). A trajectory on such a reso- 
nance chain passes from an island to the 5th next island 
and returns after r periods to the island, where it initially 
started. Considering the r-fold iterated map instead of 
the map itself, the trajectory always remains on one and 
the same island. 

The semiclassical quantization is done with respect 
to this r-fold map of period rr [25]. To each princi- 
pal quantum number m there exist r regular Floquet 
states IV^rn^) of different quantum numbers / = 0, . . . , r — 1 
with equidistant quasienergy spacing hu/r. We refer to 
these states as regular resonance states. Each of them 
has equal weights in each of the dynamically connected 
resonance islands, but with different phases. 

Similarly as in Eq. (|2T]) we derive from the semiclassical 
quasienergies the corresponding regular energies 



0.5 



^^"^ / 1 



Wm + (^)c 



(41) 



which are independent of the quantum number I. The 

(r) 

winding number Um refers to the r-fold iterated map. 

Figure [Tl^b) shows the Floquet occupations pi vs. the 
cycle- averaged energy (Ei) for the kicked rotor with 
hz = 2.35, where the phase space features in addition to 
the main regular island a 4:l-resonance around the peri- 
odic orbit of period 4, see Fig.[7fa). The entire resonance 
chain hosts 4-15 regular resonance states for h = 1/1000. 
The Floquet occupations of both the regular states of the 
central island and the chaotic states resemble those of 
Fig. [2jb). In addition, one finds a branch belonging to 
the regular resonance states. Interestingly, it has a posi- 
tive slope stemming from the fact that the cycle- averaged 
energies (Emi) of the regular resonance states Itpmi) de- 
crease with increasing quantum number m^ in contrast 
to the regular states of the central island. This is due to 
the asymmetry of the resonance tori around their respec- 
tive island center in phase space. This is another clear 
evidence, that the cycle- averaged energy does not serve 
as a suitable measure to quantify the regular occupations 
by exponential weights in analogy to the Boltzmann dis- 
tribution. 

The r regular resonance states It/Jrni) of fixed quantum 
number m have almost the same cycle- averaged energy 
{Emi)' As long as the coupling to the heat bath does 
not disturb the equivalence of the resonance islands, the 
occupations Prni of the r regular resonance states of fixed 
principal quantum number m are independent of the 
quantum number /. In Fig. [T^b) the corresponding four 
branches of the occupations Prni therefore lie almost on 
top of each other and cannot be distinguished on the scale 
of the figure. Small deviations from this degeneracy exist 
only for the outermost regular resonance states. These 
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FIG. 7. (Color online) (a) and (b) Analogous to Fig. [2] for 
the kicked rotor with ^c = 2.35 in presence of a 4:l-resonance 
chain. The insets in (b) show Husimi representations of a 
regular state from the main island, a regular resonance state, 
and a chaotic state, (c) Floquet occupations prni of the regular 
resonance states vs. regular energies ^^f, the Boltzmann- 
like prediction Eq. (|40|) with /3eff ~ 0.98/3 (red solid line) 
compared to the inverse bath temperature p (dashed line). 
The parameters are h = 1/1000 and /3 = 500. 



can be attributed to the occurrence of avoided cross- 
ings, which break the degeneracy of £mi mod (hjo/r) for 
/ = 0, . . . , r — 1 and fixed tti, as well as the degeneracy of 

Now we want to explain that the occupations Prni of the 
regular resonance states are likewise distributed as Pmi ^ 
^-(BeiiE^i , according to the exponential weights (|^Q]) with 
the effective temperature l//3eff of Eq. ([39|) . We note that 
Eq. ([29|) for the ratio of the correlation functions and 
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Eq. ()32]) for the ratio of the couphng matrix elements 
apply without restriction also to the regular resonance 
states. The assumed detailed balance relation, Eq. (|23l) . 
however, is no longer adapted to the structure of the rate 
matrix since here, in addition to the nearest-neighbor 



0.5 



rates R 



{ml){m±l,l) 



also 'internal' rates exist, i.e. rates in 



the subspace of the r equivalent regular resonance states 
/ = 0, . . . , r — 1 with fixed quantum number m. Nonethe- 
less, the total rate balance approximately decouples for 
each principal quantum number m into the r balance re- 
lations for / = 0, . . . , r — 1 



Pml 



R 



{ml){m-\-l,l) 



R 



(m+l,0(mO 



(42) 



They have the same structure as Eq. ([23|) and turn 
out to be approximately /-independent, leading to ap- 
proximately /-independent occupations Pmi • In Eq. (|^2]) 
the tiny rates R{mi){m'i') with m ^ m' and / ^ V 
are neglected and one can show that the contribution 
Er {R{mi){mi') - R{mi')imi)) vauishcs as a consequence 
of the equidistant quasienergy spacing for the r regular 
resonance states of the same m [l9[ . 

The decoupling into the r equivalent balance rela- 
tions (|^2]) finally allows us to approximate the occupa- 
tions Pml by the exponential weights e~^''^^^rni of Eq. (|4Q|) 
with the effective temperature l//3eff of Eq. (|39|) . Fig- 
ure [71(c) shows the occupations Pmi of the regular reso- 
nance states vs. the regular energies E^f. Even on the 
magnified scale of this subfigure, compared to Fig. [71(b), 
the tiny differences of the occupations Pmi with differ- 
ent quantum numbers / are not visible. The effective 
temperature l//3eff is nearly indistinguishable from the 
actual temperature 1//3, because the winding number 

v^^^jr = 0.79/4 ^ 0.2 of the resonance islands is 
small and yields a value of /3eff//3 very close to 1, com- 
pare Fig. [5l Note that it differs, although weakly, from 
/^eff ^ 0-93 of the main island. The parameter /3eff is 
the same for each of the four independent occupation 
branches Pmi • 

The phase space of a generic time-periodic system con- 
tains a hierarchy of nonlinear resonance chains and is- 
lands of all scales. If h is sufficiently small, one has Flo- 
quet states on these islands and we expect that the entire 
set of Floquet occupations becomes increasingly struc- 
tured by the branches originating from each nonlinear 
resonance. 
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FIG. 8. (Color online) (a) and (b) Analogous to Fig. [2] for 
the kicked rotor with n — 4.415 in the presence of a 2:1- 
resonance surrounded by a series of strong partial barriers. 
The insets in (b) show Husimi representations of a regular 
resonance state (vn — 3) and a beach state, (c) Floquet 
occupations prni of the regular resonance states vs. regular 
energies ^^f, the Bolt zmann- like prediction Eq. ()40p with 
/3eff ~ 0.76/3 (red solid line) compared to the inverse bath 
temperature p (dashed line). The parameters are h — 1/600 
and Q = 500. 



B. Beach states 

The transition between regular phase-space regions 
and the chaotic sea is usually not sharp, but shaped by a 
multitude of small island chains and cantori, the frac- 
tal remains of broken Kolmogorov-Arnol'd-Moser tori. 
These additional phase-space structures can strongly in- 
hibit the classical flux of trajectories toward and away 
from the regular island and, depending on the size of /i. 



can give rise to the formation of quantum beach states, 
a term introduced in Ref. [26| . These reside on the tran- 
sition layer around the regular islands and have little 
overlap with the remaining chaotic sea. Typically, beach 
states have very similar appearance and properties like 
the regular states of the adjacent island. Due to the prox- 
imity they partly even allow a quantization similar to the 
EBK-quantization rules [261.1271. 
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At K: = 4 the central island of the kicked rotor bifur- 
cates into a resonance around a stable periodic orbit of 
period 2. It is accompanied by a series of partial barri- 
ers with a reduced classical flux toward and away from 
the islands. This is indicated for hc = 4.415 in the stro- 
boscopic Poincare section of Fig. [8l^a) by the relatively 
high density of the chaotic orbit in the vicinity of the 
island. Figure Efb) shows the Floquet occupations pi 
vs. the cycle- averaged energy (Ei). The highest occupa- 
tions belong to the regular states of the resonance. The 
occupations of the beach states form a separate, nearly 
monotonous set in the transition region between the oc- 
cupations of the regular resonance states and the chaotic 
states. This is a consequence of the structure in the cou- 
pHng matrix Rij, where typically the nearest-neighbor 
rates dominate, similarly as for the regular states. 

Furthermore, the regular occupations Pmi of the regu- 
lar resonance states are shown vs. E^f in Fig. [Sfc). In 

this example the winding number u^^^q/t = 0.71/2 ^ 
0.35 in the resonance islands yields a stronger deviation 
between /3 and /3eff, with /3eff//3 ~ 0.76, than in the ex- 
ample presented in Fig. [71 
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C. Hierarchical states 

As mentioned above, in the vicinity of regular islands 
typically many partial barriers with a limited classical 
flux toward and away from the island can be found, e.g. 
in the form of cantori or based on stable and unstable 
manifolds [28] . Depending on the ratio of h to the classi- 
cal flux, partial barriers can prevent Floquet states from 
spreading over the entire chaotic domain, apart from tun- 
neling tails. If the phase-space area enclosed by the is- 
land and the partial barrier exceeds /i, these states lo- 
cally resemble chaotic states. For decreasing values of 
h they resolve and occupy the hierarchy of the classical 
phase space better and better and are therefore called 
hierarchical states [29]. The existence of these states 
does not contradict the semiclassical eigenfunction hy- 
pothesis [14], as their fraction vanishes with 0{h^) in 
the semiclassical limit. We apply an overlap criterion to 
determine, whether a Floquet state is hierarchical or not: 
it is identifled as a hierarchical state, if it is not a regu- 
lar state but comparably strongly localized, such that its 
Husimi weight fJ^dxdpH^{x^p) within a large chaotic 
phase-space area (} away from the regular island falls be- 
low 70%, compared to a state that is uniformly spread 
over the entire phase space. 

Figure M^a) shows the Poincare section and Fig. [9fb) 
shows the Floquet occupations for the kicked rotor with 
Hi = 2.5, where the fraction of hierarchical states is com- 
paratively high [29]. In Fig. [9fc), the occupations of the 
hierarchical states are emphasized. The flgure indicates 
that their occupations are distributed analogously to the 
chaotic states which explore the entire chaotic phase- 
space region. Again, the fluctuation pattern of the oc- 
cupations Pi has its origin in the randomly fluctuating 
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FIG. 9. (Color online) (a) and (b) Analogous to Fig. [2] 
for the kicked rotor with k — 2.5 in the presence of a 4:1- 
resonance surrounded by a partial barrier. The inset in (b) is 
the Husimi representation of a hierarchical state. The branch 
with positive slope belongs to the regular resonance states 
like in Fig. [Tjb). (c) Magnification of (b) with emphasized 
data points of the hierarchical states (large green crosses), 
which are determined by the overlap criterion from the shaded 
phase-space area Q in (a). The parameters are h = 1/1000 
and /3 = 500. 



rates Rij in the subspace of the hierarchical states, as is 
the case for the chaotic states. 

To conclude this section, the occupation characteristics 
of the beach states and the hierarchical states again con- 
flrm the influence of the classical phase-space structure 
not only on the spectrum and on the Floquet states, but 
eventually also on the Floquet occupations and hence on 
the asymptotic state of the system. 

Note that in the above examples. Figs. [8] and [9l either 
of the two types is predominant, but still representatives 
of the other are present. In general, hierarchical and 
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beach states coexist. For example, a few of the states 
of intermediate cycle- averaged energy {E) that are in- 
dicated in Fig. [9] as hierarchical by the above overlap 
criterion had rather to be classified as beach states or as 
states with scarring behavior, i.e. localized on hyperbolic 
fixed points or on a family of parabolic fixed points. 



VI. AVOIDED CROSSINGS 

Since the spectrum of Floquet systems is restricted to 
a finite interval < e < ficj, a multitude of avoided 
level crossings typically emerges under the variation of 
a parameter and gives rise to the hybridization of the 
involved Floquet states. In the case of an infinite dimen- 
sional Hilbert space the quasienergy spectrum is dense 
and there is no longer an adiabatic limit, i.e. any tiny 
parameter variation will hybridize infinitely many Flo- 
quet states in a complex way [13]. However, as shown 
in Ref. [12|, the asymptotic density operator p is not af- 
fected by a small avoided crossing, provided that it is 
smaller than a specific effective coupling strength to the 
heat bath. Thus, the interaction with the heat bath re- 
solves the difficulties of the dense quasienergy spectrum. 
In this section we focus on the opposite limit, where a 
single isolated avoided crossing strongly influences the 
entire set of Floquet occupations. 

Figure [10] presents a typical example of the kicked ro- 
tor. In Fig. [TOTa) the Floquet occupations pi are shown 
vs. {Ei) for two values of the kick strength near /^ = 2.9, 
very close to the parameter realization in Fig. [2l The dif- 
ference of these Ai:-values is sufficiently small, such that 
the classical phase space and almost all regular states 
vary only marginally. For two of the regular Floquet 
states, which we will denote as states a and 6, and which 
are initially identical to the semiclassical modes with the 
quantum numbers mi = 5 and m2 = 15, this is however 
not the case. Under the variation of n they undergo an 
avoided crossing dX n ^ 2.857, where they hybridize. 

The tiny /^-variation strongly affects the Floquet oc- 
cupations, and most prominently all the regular oc- 
cupations. Away from the avoided crossing {ni = 
2.85635), the regular occupations monotonously decrease 
with (^i), similar as in Fig. [2fb). When approaching 
the center of the avoided crossing (/^2 = 2.857175) the 
monotonous behaviour is locally disturbed: the states 
a and 6, as a consequence of their hybridization, have 
shifted mean energies {Eg) and (£^5) as well as modified 
occupations Pa ^ Pb- In Fig- fTOT a) the data points of the 
states a and b at hZ2 (marked as 'ac') are therefore found 
indistinguishable on top of each other. Beyond that, also 
the occupations Pm of all regular states with quantum 
numbers m from the interval [mi, 1712] change severely. 
They are close to the value of occupation pa ^ Pb of the 
hybridized states. In contrast, their mean energies (Em) 
do not change notably under the tiny a^:- variation, like 
those of the semiclassical modes mi and 7712- The rela- 
tive occupations Pm/Pm+i among the regular states with 
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FIG. 10. (Color online) Influence of an avoided cross- 
ing between the regular states ??ii = 5 and m2 = 15 on 
Floquet occupations, (a) Floquet occupations p^ vs. cycle- 
averaged energies {Ei) for a^: = Avi = 2.85 (black circles) and 
n — HI2 — 2.857175 (red dots) close to the center of the avoided 
crossing. The insets show the Husimi representations of states 
7711, 7712 at ni and of a corresponding hybridized state at H12 
('ac'). (b) Occupations pm of regular states vs. regular en- 
ergies Em^ at K — K,2 (red dots) and comparison to the an- 
alytical solution (|47|) of the rate model (solid gray line) and 
the model with Rm,m-\-i = Ro,i from 12] (dashed gray line). 
The arrow indicates the effective rate R^^ between states mi 
and 1712 according to Eq. (|44p . Note that pmi and pm2 are 
measured in the diabatic basis, in contrast to pi in (a). The 
parameters are h = 1/210 and /3 = 100. 



quantum numbers outside the range [mi, 7712] are also not 
affected. Only the absolute values of their occupations 
Pm are shifted due to the normalization J^iP* — 1- The 
latter is also the origin of a shift of the chaotic occupation 
plateau Pch- 

This example demonstrates that the presence of 
avoided crossings can change the entire character of the 
occupation distribution. To explain this impact the au- 
thors of Ref. [l2[ introduced an effective rate equation. 



= -pi^Rij^^PjRj 



(43) 
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which refers to a representation in the local diabatic basis 
of the avoided crossing, denoted by an overbar. In the 
diabatic basis the states a and b are replaced with states 
that remain invariant at the avoided crossing, i.e. the 
semiclassical modes mi and m2 in the case of an avoided 
crossing of two regular states. In Eq. (j43j) the typically 



(a) 



negligible rates Rmim2^ 
effective rate [12] 

R^ 



Rm2mi are replaced with a new 



(44) 



(hT/AY + 4^2 



which acts between the states mi and m2. The gap 
size A, i.e. the minimal quasienergy splitting \sa — 
Sbl of states a and b and the dimensionless distance 



from its center d 



{er> 



em2)/^ are characteristic 



properties of the avoided crossing. Unlike the rates 
Rij^ which are nearly constant in the vicinity of the 
avoided crossing, the additional rate R^^ changes dra- 
matically. The composite rate T = ^ • {Rmij + ^m2j) ~ 
{27T/h) Y^j^ Armmi {K)Am2m2 {K)g{-Khw) plays the role 
of an effective coupling strength and the characteris- 
tic parameter hT/A determines, whether R^"^ can be- 
come dominant at the center of the avoided crossing 
d = 0. In the examples of Figs. [10] and [11] the condi- 
tion hT/A <C 1 is fulfilled and the rate i?^^ consequently 
dominates around d ~ with respect to other rates in 
Eq. (j33|). Note that this requires the system-bath cou- 
pling strength to be sufficiently small, a limit which has 
been already presumed in Eq. (|6|). The dominance of 
^ac ^g responsible for the local disruption of the formerly 
exponential behavior of the regular occupations: as R^^ 
exceeds all other rates, it induces occupation equality of 
the states m^i and 7712- The relative occupations of the 
regular states below mi and above m.2 still follow the 
approximate detailed balance, Eq. (|23]) . with rates act- 
ing predominantly between nearest neighbors. Finally, 
the other regular states, m^i < mt < 7712, although also 
still having dominant nearest-neighbor rates, no longer 
have exponentially scaling occupations, since the addi- 
tional rate R^'^ between mti and r?i2 breaks the detailed 
balance. These conditions explain the observed signature 
of the avoided crossing in the occupation characteristics 
and are substantiated by a simplified rate model intro- 
duced in the following section. 

A second example, with an avoided crossing between 
the regular state mi = 10 and a chaotic state, is pre- 
sented in Fig. [TlT a). Here, the occupation of the state 
mil is forced down to the chaotic occupation level in the 
vicinity of the avoided crossing (red dots). Moreover, 
the entire subset of occupations Prn with m > mi is dis- 
turbed from the original exponential scaling (black cir- 
cles). With a similar reasoning as above, this behavior 
is explained by the additional rate R^'^ and a simplified 
rate model, see below. 

These consequences of avoided crossings also explain 
why the transition between the occupations of regular 
and chaotic states appears so 'smoothly', see e.g. Figs. [1] 
and [2] As the outer regular states or the beach states 
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FIG. 11. (Color online) Influence of an avoided crossing 
between the regular state mi = 10 and a chaotic state on 
the Floquet occupations, analogous to Fig. [TO] for k = ki = 
2.856400 (black circles) and k = K2 = 2.856897 (red dots) 
close to the center of the avoided crossing. Note in (b), that we 
have simulated the avoided crossing of ??ii with a chaotic state 
in the rate model, Eq. ()47p . by an avoided crossing between 
mi and m = 22, the last regular state. The parameters are 
h = 1/210 and 13 - 100. 



most likely undergo avoided crossings with chaotic states, 
the occupation probabilities normally feature a transition 
to the plateau pch of the chaotic occupations. 



We emphasize that the observed implication of an 
avoided level crossing is a remarkable effect bound to 
the non-equilibrium character of the driven system, with 
possible applications, e.g. bath-induced switching in a 
driven double well potential [ISj. In a time-independent 
system on the contrary, an avoided crossing entails only 
a local shift in the occupations of the two involved states 
and leaves the Boltzmann distribution of the entire set 
of occupations unchanged. 
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A. Simplified rate model 

To account for the local modifications of the occupa- 
tions Pm^ which are entailed by an avoided crossing, we 
apply a simplified, analytically solvable model for the 
set of rate equations (j43j). It is based on the nearest- 
neighbor rate model in Ref. [12|- It is restricted to the 
chain of regular states and, at first in the absence of 
an avoided crossing, assumes that each of them is cou- 
pled only to its directly neighboring states. Under this 
circumstance the detailed balance is fulfilled. The pri- 
mary parameter of the model is the rate ratio of the 
first two states, g := i^o,i/^i,o- We approximate the 
rate ratio Rm,m-\-i/Rm-\-i,m as independent of m, i.e. 
Rm,m-\-i/ Rm-\-i,m = 9 for all m. According to Eqs. (|38|) 
and ([23|) the approximation is especially suited for regu- 
lar islands with a constant winding number v. 

The rates themselves however differ from state to state. 
In contrast to Ref. [12], where Rm,m-\-i/Ro,i is presumed 
as constant, we make use of the approximation 



il-9"''-"'')[rmJil+r^,)] +9" 



approaches 1. 



R 



m,m+l 



^0,: 



m - 



1. 



(45) 



This relation holds exactly for the states of a harmonic 
oscillator with a linear coupling to the heat bath. The 
coupling matrix elements fulfill Am,m-\-i = V^ + 1 ^0,1, 
a property that is easily shown with the help of the as- 
sociated algebra of ladder operators. Analogously one 
can prove relation (|45]) for the regular states of a time- 
periodic system if the regular tori are elliptic, coincid- 
ing with those of a harmonic oscillator, and assume m- 
independent winding numbers. Beyond the application 
to elliptic islands of constant winding number we pre- 
sume Eq. (|45]) for arbitrary islands, which seems to be in 
general a good approximation. 

To account for an avoided crossing between the states 
mi and 7712 the model introduces the additional rates 
Rmim2 = Rm2mi = R^^ from Eq. (|44|) . By the action 

of i?^^ the fiux Fjn^rn-l •= PmRm,m-l — Pm-lRm-l,m 

becomes non-zero for all mi < m < m2. The rate equa- 
tions, Eq. (|43|) . translate to the fiux equations 



= Fin 



(46) 



— Fm-^i^rn — ^m,m-l m^0,mi,m2 

[Pmi Pm2 ) -^ — -^mi+l,mi -^mi,mi — 1 



\Pmi Pm2 ) -^ — -^7712 + 1,' 



F, 



m2,m2 — 1 



with the solution Fm,m-i = F := {pmi — Pm2) ^^"^ fo^ 
mi < m < 7712 and Fm,m-i = otherwise. Eventually, 
the occupations assume the values 



P 



Po . 

Prm 
Pm2 9 



(1-^^2-mi) 
,m—m2 



l+r„ 



mi<7n<7Ti2 
m2<m 

(47) 

with Tm := i^^^ELT' (^"'""'^"V^m,+fc,m,+fc-l). For 

^ac ^ ^^^ |-]^g parameter r^2 diverges and Pm2/Pmi = 



Note that the model solution (|T7|) relies on the 
nearest-neighbor coupling and on the m-independence of 
^m,m+i/^m+i,m. The ratio Rm,m+i/Ro,i of Eq. (gSl) 

leads to Tm = (i^^7^i,o) Er=T' [^^-^^- V(^i + k)] . 

In Figs. [TOlf b) and \TT\ h) we apply the simplified rate 
model to regular states for a kick strength close to the 
center of the avoided crossing, and compare its solution, 
Eq. (|47|) . to the occupations Pm from the solution of the 
rate equations (|6]) . The regular energies E^^^ are not well 
defined for the hybridizing states of the avoided crossing, 
but are well defined for the respective diabatic states mi 
and 77i2- That is why the occupations in these subfigures 
are represented in the diabatic basis by means of the 
orthogonal transformation p^^ = a'^Pa^^'^Pb andprn2 = 
P'^Pa^a^Ph with Qf^ = (1 + d/y/l^d?) /2 and a^ + Z^^ = 
1. The comparison indicates that the rate model based 
on assumption (|45]) reproduces the local disturbance of 
the exponential scaling for the states mi < m < m2 very 
accurately. In contrast, the simpler assumption [12| of m- 
independent rates, Rm,m-\-i/Ro,i = 1^ does not reproduce 
the 777-dependence of the occupations between m^i and m.2 
(dashed line). 

For the example in Fig. [TlTb) the rate model seems at 
first sight not applicable, since it does not account for 
chaotic states. We therefore apply the model as if the 
avoided crossing were between the regular state mi = 
10 and the outermost regular state, m = 22, and still 
obtain a good agreement between the observed pm and 
the model solutions. This is possible, as the occupation 
of the regular state tti = 22 differs only weakly from the 
plateau of chaotic occupations. 



VII. SUMMARY 

A core question of statistical mechanics is the char- 
acterization of the asymptotic state approached by a 
quantum system when it interacts with a thermal reser- 
voir. In the familiar equilibrium thermodynamics of 
time-independent systems in the weak-coupling limit it 
is answered by the canonical distribution, where the 
eigenstates of the isolated quantum system are occu- 



pied with the statistical weights pi 



-I^Ei 



In a time- 



periodic quantum system, where an external field per- 
manently pumps energy into the system and prevents 
its relaxation to equilibrium, this is in general an intri- 
cate question, which cannot be answered by deduction 
from the time- independent case. Here, the asymptotic 
state under a weak coupling to the thermal reservoir 
becomes time-periodic and is best characterized by the 
time-independent occupations pi of the Floquet states. 
We demonstrate, that the Floquet occupations can be 
classified according to the semiclassical character of the 
Floquet states. The occupations of the chaotic Floquet 
states fiuctuate weakly around a mean value Pch [lO]- 
The regular Floquet states on the contrary acquire prob- 
abilities that are roughly exponentially distributed. The 
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validity of this observation is also confirmed by the oc- 
cupation characteristics of other types of Floquet states, 
which still reflects their regular or chaotic nature: beach 
states, which are very similar to the regular states and 
situated close, but outside the regular island, they form a 
correlated set of occupations, which is qualitatively com- 
parable to the regular occupations. In contrast, the occu- 
pations of hierarchical states, which have the properties 
of chaotic states, but live in a restricted region of the 
chaotic phase space, are distributed analogously to the 
chaotic states. 

In contrast to previous studies of a driven particle 
in a box [10], where the regular states carry occupa- 
tions close to the Boltzmann weights, we observe that 
in general the regular occupations can considerably de- 
viate from the Boltzmann result. This observation is 
possible as we focus on time-periodic systems where the 
classical phase space and the Floquet states are strongly 
perturbed compared to the originally time-independent 
system. In kicked systems the occupations of the reg- 
ular states can be well approximated by a function 
Fivrn-, P^) depending on the classical winding numbers 
Vrn of the regular tori, and the parameters temperature 
1//3 and driving frequency uo. For a constant or suffi- 
ciently moderately varying winding number within the 
regular island the distribution of the regular occupations 
can even be described by weights of the Boltzmann type, 
Pm ^ e~^^ff^^ , depending on the regular energies £^^^. 
The effective temperature l//3eff is evaluated as a func- 
tion of the winding number in the regular island. As the 
driven system is not in an equilibrium state, the proper 
definition of a non-equilibrium temperature is a subtle 



problem [30|. The critical question for future investi- 
gations is, whether the quantity /3eff is accessible by a 
measurement. 

A situation, where purely classical information is no 
longer sufficient to account for the observed occupations, 
is present at avoided crossings, which are ubiquitous in 
the quasienergy spectra of Floquet systems. Avoided 
crossings involving a regular state give rise to strong 
changes in the set of Floquet occupations. We give an in- 
tuitive explanation of this effect, based on the additional 
rate R^^ from Ref. [12], and introduce a simplified rate 
model whose analytical solutions describe the numerical 
data accurately. 

In conclusion, with the presented characterizations of 
the Floquet occupations we demonstrate that ubiquitous 
signatures of the classical dynamics are reflected in the 
asymptotic density matrix of the open quantum system. 
In this way it is feasible to draw an intuitive picture of 
the asymptotic state, shedding light on the statistical 
mechanics of time-periodic quantum systems. 
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